Systematic study of low-lying El strength using 
the time-dependent mean field theory 

S. Ebata*' + , T. Nakatsukasa^** and T. Inakura 1 ^ 

* Center for Nuclear Study, Univ. of Tokyo, Bunkyo-ku, 113-0033, Japan 
1 Theoretical Nuclear Physics Laboratory, RIKEN Nishina Center, Wako-shi, 351-0198, Japan 
** Center for Computational Sciences, Univ. of Tsukuba, Tsukuba-shi, 305-8571, Japan 

Abstract. We carry out systematic investigation of electric dipole (£1) mode from light to 
heavy nuclei, using a new time-dependent mean field theory: the Canonical-basis Time-Dependent 
Hartree-Fock-Bogoliubov (Cb-TDHFB) theory. The Cb-TDHFB in the three-dimensional coordi- 
nate space representation can deal with pairing correlation and any kind of deformation in the time- 
dependent framework. We report the neutron-number dependence of the low-energy E 1 mode for 
light (A < 40) and heavy isotopes (A > 100) around N = 82. 
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INTRODUCTION 

Properties of the low-energy electric dipole (El) excited states of neutron-rich nuclei 
are important to understand nucleosynthesis on the r-process path [1]. These states are 
often called Pygmy dipole resonance (PDR). Due to the recent progress of radioactive 
facilities, we become able to generate and measure the neutron-rich nuclei, however 
it is still difficult to study nuclei on the r-process path. Thus, we strongly demand 
a method (model) that is able to analyze and predict the excitation modes of exotic 
nuclei in very neutron-rich and heavy nuclear region, such as nuclei on the r-process 
path. In order to study dynamical properties of heavy unstable nuclei, the method 
should take into account the effects of deformation and pairing correlation. The Time- 
Dependent Hartree-Fock-Bogoliubov (TDHFB) or HFB plus quasi-particle random- 
phase approximation (HFB+QRPA) are such candidates [2], which are applicable to 
a wide range of nuclei from light- to heavy-mass regions. However, they require a 
significant effort for coding the program as well as large computational resources. The 
TDHFB calculations with full three-dimensional (3D) dynamics is currently still in a 
preliminary stage [3,4]. We propose a feasible approach, AgCanonical-basis TDHFB Ah 
(Cb-TDHFB) in the 3D coordinate- space representation [5]. The numerical cost of Cb- 
TDHFB is much smaller than that of TDHFB, then we can study excited states of various 
isotopes systematically. 

In this present work, we investigate the low-lying E 1 modes of even-even neutron- 
rich nuclei using the Cb-TDHFB in the linear regime. We report that the low-lying El 
strengths of light and heavy neutron-rich isotopes have characteristic neutron-number 
dependence. 



FORMULATION 



The Cb-TDHFB equations [5] are derived from the TDHFB equations with a simple 
approximation for the pairing interaction which is using only diagonal part of pairing 
functional A/£ = A/ 8j k analogous to the BCS approximation. In present work, we choose 
a very schematic pairing functional as £p a i r = —Y.i,m>oGim K l K m with fQ = w/v/. u\ 
and v/ correspond to the time-dependent BCS factors for the canonical pair of states, 
<j>l(r, a,t) and <^-(r, a,?). The Cb-TDHFB equations for the canonical single-particle 
states, the occupation probabilities p t (t) = |v/(?)| 2 , and the pair probabilities K t (t) are 
written as 

{ffi<j>i(r,o,t) = (h[p(t)} -£/(O)0z(r,(7,O, in(j>i(r,a,t) = {h[p(t)} - e^))0f( r ><M)> 
ifip l (t) = K l {t)A^t)-A l (t)Kf(t), 
iHki{t) = (£i(t)+e T (t))Ki(t) +A/(0(2p,(0 - 1), 

where, /?[p(?)] is the HF Hamiltonian. A/(?) and £/(?) are defined as A/(?) = 
Hm>oGim K m(t) and £/(?) = respectively. The Cb-TDHFB equa- 

tions conserve the orthonormal property of canonical basis, the expectation value of 
particle number and total energy. They are identical to TDHFB. We use the Skyrme 
functional of the SkM* parameter set for /^-channel and G/ m = G/(e / °)/(e^), here 
/(e) is a cutoff function, adopted from Ref.[6] for pp,hh-channe\. 

In the liner-response calculation, we add an external field which is weak and instan- 
taneous in time, V ex t(r,t) = —kF(r)8(t). F is one-body operator. In this work, we adopt 
the £1 operator, F = (Ne/A)r^ — (Ze/A)f^ n \ where r = (x,y,z) and (p), (n) mean the 
operation on protons and neutrons, respectively. If the parameter k is very small, the 
fluctuation of the nuclear density is linearized automatically, and this linear-response 
calculation is equivalent to the QRPA (or RPA). We obtain the strength function S(F;E) 
through the Fourier transformation of the expectation value of the one-body operator F. 
Details can be found in Refs.[5, 7]. 



RESULTS 



We discuss a systematic trend of the low-lying El strength. To quantify PDR, the 
following ratio of the low-lying E 1 strength is used: 

nnm S E >E*S(EV,E)JE 
mi fExS(El;E)dE 1 h V 

where E c is a cut-off energy. We take E c =10 MeV in the present calculation. 

Figure 1 shows the m ^/ > f° r O, ^ e and Mg isotopes as a function of neutron number. 
We see that the nuclei with /V =8 - 14 have the ratio of low-lying El strength less 
than 1.0%, and a sudden jump of the ratio at /V = 14 — > 16 on each isotopic chain. The 
neutron number N = 16 corresponds to the occupation of the s\ji orbit. These results 
of light neutron-rich nuclei suggest that the position of the neutron Fermi level plays an 
important role for the emergence of low-lying El strength [8]. 



Figure 2 shows the results for heavy neutron-rich isotopes which are Z = 38 — 54 
with N = 76 — 90. Again we can see that a sudden jump of the ratio at N = 82 — > 84 
on each isotopic chain. These results show a neutron shell effect on the low-lying El 
strength. However, there are constant but significant PDR ratios for N = 16 — 82. Al- 
though 13 0-i36x £ are stable nuclei, they have 1.5% of low-lying El strength. This ratio 
systematically increases as the proton number decreases. These features suggest that the 
PDR in heavy nuclei has a different structure from one of light nuclei. Currently, we are 
investigating the origin of the difference, and also studying even heavier neutron-rich 
nuclei around N = 126. 




Neutron Number 

FIGURE 1. Neutron-number dependence of m\ (E c ) /m\ defined in Eq.(l) for O, Ne and Mg isotopes. 
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